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Abstract 

Motivated by recent investigations of transport properties of strongly correlated 
Id models and thermal conductivity measurements of quasi Id magnetic systems we 
present results for the integrable spin-1/2 XXZ chain. The thermal conductivity k(uj) 
of this model has Re«(w) = k8(u>), i.e. it is infinite for zero frequency to. The weight k 
of the delta peak is calculated exactly by a lattice path integral formulation. Numerical 
results for wide ranges of temperature and anisotropy are presented. The low and high 
temperature limits are studied analytically. 
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1 Introduction 

Transport properties of strongly correlated quantum systems have recently attracted strong 
theoretical and experimental interest. For the XXZ chain as a quantum spin model (or 
alternatively as a lattice gas model) the spin transport (electrical transport) was investi- 
gated analytically [|J by use of a method suggested in Q and numerically ||,|4) with so far 
inconclusive results about the finite temperature Drude weight. In the thermal con- 
ductivity of the quasi one-dimensional magnetic system (Sr,Ca)i4Cu2404i was investigated 
showing anomalous transport properties along the chain directions. In (?J it was argued 
that the large thermal conductivity can not be explained in terms of phonons. 

Here we want to use an approach as microscopic as possible and apply Kubo's theory to 
the strongly correlated spin-1/2 XXZ chain. This approach closely follows ||. We explic- 
itly avoid the notion of particles and Boltzmann equations that take account of scattering 
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as is commonly done in Fermi liquid theory. In strongly correlated low-dimensional quan- 
tum systems there is neither a particle picture with finite quasi-particle weight nor do we 
restrict ourselves to low temperatures in comparison to some reference ("Fermi") energy. 
Rather on the opposite, we are interested in temperatures comparable to the exchange 
energy of neighboring spins. 

The Kubo formulas ||[l(J are obtained within linear response theory and yield the 
(thermal) conductivity k relating the (thermal) current J7e to the (temperature) gradient 
VT 

Je = kVT, (1.1) 

where 



/ dte- iujt / dT{J E (-t-iT)J E ), (1.2) 
Jo Jo 



and (3 is the reciprocal temperature 1/(&;bT)[]. The current-current correlation function is 
to be evaluated in thermal equilibrium and poses the main problem of our work. However, 
as already pointed out in [8, ll-15fl the total thermal current operator J7e commutes with 
the Hamiltonian TL of the XX Z chain. Hence we find 

1 

\{uj — ie) 



k(u) = -. 1 . , P 2 (Ji), (e -» 0+), (1.3) 



with Re = k5(uj) where 

K = ir(3 2 (Ji). (1.4) 

As a consequence, the thermal conductivity at zero frequency is infinite! The property 
of a non-decaying thermal current is a feature of Id systems in particular of integrable 
systems with many integrals of motion (for a discussion see fl23"|| ). We expect that additional 
interactions, especially residual interchain couplings in quasi one-dimensional materials will 
lead to a broadening of the delta peak, however with same weight. Instead of speculations 
about the possible scenarios of such a broadening we want to present results for the weight 
k as a function of temperature and anisotropy. 

The paper is organized in the following way. In Sect. 2 we discuss the integr ability of the 
XX Z chain and identify the thermal current as one of the integrals of motion. In Sect. 3 we 
present our method of calculation of thermodynamical quantities and discuss the results. 
In Sect. 4 we treat analytically the low and high temperature limits and give a summary in 
Sect. 5. 

2 Thermal current and conserved quantities 

Here we begin the microscopic approach to the heat conductivity of the spin- 1/2 quantum 
chain with a construction of the conserved currents. The Hamiltonian of the XX Z model 



1 In this paper we set fes = 1. 
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on a periodic lattice of size L is 



L 



fc=i 



hkk+i = J ( o-fc <r k+1 + o-fc +1 o- fc + iro-fccrf+i ) > (2- 1 ) 



where cr^ = (cr? ± ia"|)/2 and cr^, cr|, denote the usual Pauli matrices acting on the kth 
space. 

In this paper we restrict ourselves to the critical regime — 1 < A < 1 where the system 
displays algebraically decaying correlation functions at zero temperature. The anisotropy 
parameter A is conveniently parameterized by 

A = cos 7, < 7 < 7T. (2.2) 

We will give our analytical results for the entire regime < 7 < 7r (— 1 < A < 1), numerical 
results will be given for the repulsive range < 7 < tt/2 (0 < A < 1). 

Our first goal is the determination of the thermal current operator j E . To this end 
we impose the continuity equation relating the time derivative of the local Hamiltonian 
(interaction) and the divergence of the current: h = — div j E . The time evolution of the 
local Hamiltonian in ( [2.1D is obtained from the commutator with the total Hamiltonian 
and the divergence of the local current on the lattice is given by a difference expression 

= m h kk+1 ( t )\ = -{jf +1 (t)-jf(t)}. (2.3) 

Apparently the last equation is satisfied with a local thermal current operator j E defined 
by 

Jk = i[hk-ik,hkk+i]- (2-4) 
In fact the total thermal current J~e = Y^k=iik 

L 

J E = -iJ 2 J2{^ti a k+1 ~ a k+X a k-i) ~ A K-i + frfc+aX^fc " °t+i°k)}' ( 2 - 5 ) 
fc=i 

commutes with the Hamiltonian, [H, ^7e] = 0, as it is closely connected with a nontrivial 
conserved quantity derived from the underlying integrability of the model. To see this, we 
introduce the transfer matrix constructed from the ii-matrix R £ End(F (8> V) where V 
denotes a two dimensional irreducible module over the quantum algebra U q (sl(2)). The 
nonzero 6 elements of the i?-matrix with the spectral parameter v are given by 

R\\(v) = RjUv)= [ ^, R\t(v) = Rl\(v) = R%(v) = R£(v) = l, (2.6) 
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where [v] denotes [v] = sm{^v/2)/ sin (7/2) and the indices of the above relations can be 
interpreted as 

|i) = IT>, |2) = U>, 

R{v)\a)®\(3)=Y,\l)®\5)&i{v)- (2-7) 

7,5 

The i?-matrix satisfies the Yang-Baxter equation (YBE) 

R 12 (u - v)R 13 (u)R 23 (v) = R 23 (v)R 13 (u)R 12 (u - v). (2.8) 

Due to the YBE the transfer matrix T(v) = Y\ k Rkk+i( v ) 1S commutative with respect to 
different spectral parameter v and v' , i.e 

[T(v),T(v')]=0. (2.9) 

As is well known, \nT(v) is the generating function for the conserved quantities 

J (n) = (-P\ n ^T(v) . (2.10) 

In particular the Hamiltonian ( |2.1| ) and the thermal current ( |2.5| ) are expressed in terms 
of and J^ 2 \ respectively. Explicitly they read 



n = 2 l^ljm - —A, 
7 2 

J E = i(^^) 2 J^ + U*L. (2.11) 
Due to the commutativity (|2.9| ), every operator J'^- n > commutes with the Hamiltonian 7i. 

3 Thermal conductivity 

Our goal is the calculation of the second moment of the thermal current. Quite generally the 
expectation values of conserved quantities may be calculated by use of a suitable generating 
function. As such we want to define a modified partition function 

Z = Tr exp(-(3H + \J E ), (3.1) 
from which we find the expectation values by derivatives with respect to A at A = 



a=o V d\ 



a x2 



m = 0, — InZ = (Ji) - (J E ) 2 = (J4), (3.2) 



A=0 



where we have used that the expectation value of the thermal current in thermodynamical 
equilibrium is zero. 
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Instead of Z we will find it slightly more convenient to work with a partition function 

Z = Tr eM-^iJ {1) ~ KJ {n) )- (3.3) 
With view to (2.11) we choose 



Ai = /3 



2Jshi7 

7 



A 2 = — iA 



2Jsin7 

7 



(3.4) 



and obtain the desired expectation values from Z 



d 



A=0 



(3.5) 



We can deal with Z rather easily. Consider the trace of a product of N row-to-row 
transfer matrices T{uj) with some spectral parameters Uj close to zero, but still to be 
specified, and the Nth power of the inverse of T(0) 

Z N = Tr [T(m) T(u N ) • T(0)~ N ) 

= Tr exp Yji lnT ( u j) ~ lnT (0)] 



Now it is a standard exercise in arithmetic to devise a sequence of N numbers ui,...,un 



(actually Uj = u^) such that 



d 



dv' 



v=0 



v=0 



(3.7) 



We only need the existence of such a sequence of numbers, the precise values are of no 
importance. In the limit N — > oo we note 



lim Zjv = Z. 

N^>oo 



(3. 



We can proceed along the established path of the quantum transfer matrix (QTM) 
formalism |TB-p^] and derive the partition function Zn in the thermodynamic limit L — > oo 



lim Z 



1/L 
N 



A. 



(3-9) 



where A is the largest eigenvalue of the QTM. The integral expression for A reads 

/oo 1 
K{v) HK{v)K(v)]dv, K(v) = 4cQgh ^ , (3.10) 
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with some function e(v) given in |21|. In the limit N — > oo the first term on the r.h.s. of 
the last equation turns into 







N^oo 



v=0 



e(v 



v=0 



(3.11) 



a rather irrelevant term as it is linear in Ai and A n , and therefore the second derivatives 
with respect to Ai and A n vanish. The functions 21(d) and 21(f) are determined from the 
following set of non-linear integral equations (NLIEs): 



lna(v) = ^2[e (v - iuj) - e (0)] + k * ln2l(u) - k * ln2l(t> + 2i), 
j 

hxa(v) = ^2[s (v - iiij) - e (0)] + k * In 21(f) - « * ln2t(f - 2i), 

3 

a(v) = l + o(u), §f(«) = l + o(u). 



(3.12) 

with a function £o( u ) given in terms of hyperbolic functions pi] . The symbol * denotes 
the convolution / * = /(« — v')g(v)dv and the function is defined by 



k(v) 



1 

2^ 



sinh 



2 fc 



oo 2 cosh fc sinh f ? — 1 ) 



(3.13) 



Again, the summations in ( |3.12| ) can be simplified in the limit N — > oo 



lim y^[eoi 



f — ltlj 



eo(«)] 



) -<><» -A, 



:ei(u) 



-l^J ^o(f) 



■■s„(v) 



(3.14) 



where the first function can be found in [21] and is simply 



s x {v) = 2nK(v) 



TT 



2 cosh ' 



and hence the second function is 

e n (v) 



,d_ 

dv 



n-l 



ei(v). 



(3.15) 



(3.16) 



We like to note that the structure of the driving term ( 3.14 ) appearing in the NLIE fljDg ) 
reflects the structure of the generalized Hamiltonian in the exponent on the r.h.s. of (|3,3j). 
We could have given an alternative derivation of the NLIE along the lines of the ther- 
modynamic Bethe Ansatz (TBA). In such an approach the driving term is typically the 
one-particle energy corresponding to the generalized Hamiltonian. Hence it has contribu- 
tions due to the first as well as the nth logarithmic derivative of the row-to-row transfer 
matrix, i.e. the terms e\ and e n . 
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In FigjT] (a) we show H(T) for various anisotropy parameters 7. Note that k(T) has 
linear T dependence at low temperatures. In the free fermion case (7 = it/ 2) we observe 
that the ratio of the thermal conductivity k to the electrical conductivity a obeys the 
Wiedemann- Franz law (see the next section). There is a finite temperature maximum at 
roughly half of the temperature of the maximum in the specific heat c(T) data. At high 
temperatures ii(T) behaves like 1/T 2 . This and the low temperature asymptotics will be 
studied analytically in the next section. 

In general the data of h{T) show a much stronger variation with the anisotropy 7 than 
the specific heat c(T) data. In Figjl] (b) we show the ratio k(T)/c(T) which have finite low 
and high temperature limits and strongly depend on the anisotropy parameter. 



4 Low and high temperature limits 

In certain limits we can analytically evaluate the NLIE. For high temperatures the NLIE 
linearize and hence can be solved. At low temperatures the NLIE remain non-linear, but 
the driving terms simplify. Then the symmetry of the integration kernel allows for an 
evaluation of the physically interesting quantities by avoiding the explicit solution of the 
NLIE. 

4.1 Low temperature limit 

We consider the low-temperature behavior of the current-current correlation function of 
the thermal current by making use of the dilogarithm trick. The functions a(v) and 
(i(v) in the NLIEs (pit ) exhibit a crossover behavior 



a(v) < 1, a(v) < 1 for -K- < v < K+, 

a(v) ~ 1, a(v) ~ 1 for v < -K-, K + < v, (4.1) 

where 

K ± = 2 - In { (b ± I . (4.2) 

We introduce the scaling functions 

a±(v) = a (±^v ±/C±^ , a±(v) = a(±^v±K±\ , (4.3) 

and similarly the functions A± (^4±) for 21 (21). From the NLIEs ( |3.12j ), one sees the scaling 
functions satisfy the following "scaled" equations in the limit (3 — > 00: 

lna±(v) = —e~ v + K\ * lnA±(v) — K2± * lnA±(v), 

lna±(v) = —e~ v + K\ * ]nA±(v) — K2 T * lnA±(v), 

2 (2v\ 2 (2v \ , s 

«!(«) = -« — , k 2 ±(v) = -k — ±2i . (4.4) 



IT \ IT J 7T \ IT 
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In this limit, the integrals of the functions 21, 21 in eq. ( |3.10| ) can be written as 

/oo 
K{v) Yn[%(v)$L(v)]dv 
-oo 

= 9 2 t ■ TW~L XT ■ S I" e-"(ln^ + ( W ) +lnA + ( W ))d« 
27r" ! Jsm7(p7 + A7rJsm7j 

+ o — ,f x T ■ \ r e- v {lnA„(v) + lnA.{v))dv. (4.5) 



27T 2 Jsin7(/?7 — A7rJsin7) 

The right hand side of the above equation is evaluated as follows. (A) After taking the 
derivative of the first and second equation in (|4,4|), we multiply them by lnA±(v) and 
In ^4-|- (u), respectively and take the summation of them with respect to each side. (B) Next 
we multiply (4.4) by [ln^4.±(u)]' and [In A±(v)]', respectively and take the summation. Next 
we subtract the resultant equation of (B) from the one of (A). After integrating over v, we 
obtain 

/oo 
e- v \kiA±(v)+hiA±(v)]dv, (4.6) 
-oo 



where 



f°° ( d — d 

D± = / (]nA±(v)—\na±(v)+lnA±(v)—lna±(v) 
J-oo v d v dx 

— In a± M — lnA±(v) — lna±(v) — In A±(v))dv 
dv dx J 

a±(oo) /l n (l +a ) ma \ fa±{oa) /U 1+ j) fog 

1 da + / — — - - = da. (4.7) 



/a±(-oo) V a 1 + a/ Ja±(-oo) V a l + a 

The quantities -D± can be expressed in terms of Roger's dilogarithm C(v 



l + a±(oo)/ \l + a±(oo 



1 + a±(— oo)/ \ 1 + a±(— oo) 

1 r /ln(l -x) Lax \ 

Using the asymptotic value of the scaling functions a±(±oo) = a±(±oo) = 1 and substi- 
tuting eq. ( fO| ) for eq. (4J5), we arrive at 

lnA = - [ - | - V (4.9) 

12/3Jsin7 \/?7 + A7rJsin7 /?7 — A7rJsin7/ 

Here we have used the identity C(v) + C{1 — v) = 7r 2 /6. Hence the low-temperature 
asymptotics of the current-current correlation function is evaluated to 

«2> = ^W?)- 
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From this result we see that k{T) is linear in T at low T 

_2 



7T 



k{T) ~ — i/T, 



(4.11; 



with u = j7rsin7/7 the velocity of the elementary excitations. From the low temperature 
behavior of the specific heat c(T) = {tt/2>v)T we find 



c(T) 



(4.12) 



Finally we want to compare the thermal conductivity k and the electrical conductivity a 
which in the low temperature limit can be described by the Drude weight D c 



Rea(u) = 2ttD c 5(lo), D c 



4(7T - 7) 



This yields 



-~|tt (tt-7)T, (T^O), 
a o 



(4.13) 



(4.14) 



which in the free fermion case (7 = 7r/2) gives the Wiedemann- Franz law. 
4.2 High temperature limit 

In the high temperature limit {(3 — > 0), the auxiliary functions satisfy the following integral 
equations linear in (9/9A) 2 ln2l 

A) 2 ln2l(u) = (^m 2 t ( ,)) 2 + 1.* (A) 2 ln2l(u) _ ^* ^\ n ^ v + 2 i), 



d . = 



9A 



ln2l(u 



1 / 9 



~ 1 / 3 \ ^ 

m2t(u) ( — ] ln2l(v-2i). 
2 \oX I 



(4.15) 

Here we have used the high temperature asymptotics a(v),a(v) ~ 1. By use of the dressed 
function formalism, we obtain the identity 



A=0 



2T / (l^) \n*icmc) 



7 



8ttJ sin 7 



9A 



<9a(u) / da(v 



dv 



A=0 



9/3 



OA 



dv. 



A=0 



(4.16) 



The integrand in the above equation is found analytically 



da(v) 



d(3 
da(v) 



J sin 2 7 



2 sinh ^(t> + i) \ sinh ? (f + 3i) sinh^(?;— f 



9 / da(v) 



=0 <9v \ d/3 



(4.17) 



A=0 
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By use of these explicit expressions we obtain the high temperature limit of k(T) 

k{T) S -^) (3 2 + 0(/? 3 ). (4.18) 



4 \ sin 7 

This result generalizes the known results for the special cases 7 = tt/2 (free fermion model) 



in 1 25 1 and 7 = (isotropic Heisenberg chain) in |j26[ . 



5 Summary and discussion 

In this work we have presented a method for the calculation of the thermal conductivity 
k(lu) of the integrable spin-1/2 XX Z chain. The investigation of this system is drasti- 
cally simplified in comparison to other systems as here the thermal current is a conserved 
quantity. This led to Reft(w) = k5(lv) where the weight k was calculated for arbitrary tem- 
perature and various anisotropy parameters in the repulsive, critical regime of the XX Z 
model. 

We like to comment on possible generalizations of our investigation. First of all, the 
computation of K,(q,u>), i.e. the thermal conductivity for a thermal current operator Je(q) 
with non-zero momentum q, would be desirable. Unfortunately, it is only the case q = 
that allows for an analytical approach. On the other hand, for small values of q a behavior 
of n(q,Lo) similar to the q = case is expected with the 5(u>) factor to be replaced by 
5(uj — vq) where v is the velocity of the elementary excitations. Finally, the investigation of 
the thermal conductivity of the general XX Z chain with arbitrary anisotropies in the easy 
plane (gapless) and easy axis (gapped) regimes should be possible. This will be reported 
elsewhere. 
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Figure 1: (a) Illustration of numerical results for the thermal conductivity h as a function of 
temperature T for various anisotropy parameters 7 = 0, tt/Q, n/5, 7r/4, 7r/3, tt/2. (b) Depic- 
tion of the ratio of thermal conductivity and specific heat k/c as a function of temperature. 
The analytic result ( 4.12j ) in the low temperature limit is also depicted by the symbol *. 
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